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39. Current Evolution of Meteoroids 



J. S. Dohnanyi 
Bellcomm, Inc. 
Washington, D.C. 


The observed mass distribution of meteoroids at 1 A U from the Sums briefly reviewed 
in a survey that ranges over the bulk of the mass spectrum from micrometeoroids to 
meteorite parent objects. The evolution of meteoroids under the influence of collisions, 
planetary perturbations, the Poynting-Robertson effect and radiation pressure is then 
discussed. 

Most micrometeoroids are expelled from the solar system by radiation pressure shortly 
after their production as secondary ejecta during impact by larger objects or as dust 
ejected by comets. Particles that survive will eventually be swept out by the Poynting- 
Robertson effect. 

Meteoroids in the radio and photographic ranges are destroyed in collisions faster 
than they can be replaced by the production of secondary fragments during collisions 
bet ween larger objects . The source of new particles needed to maintain the population of 
these meteoroids in a stationary distribution may be material expelled by comets. 

The survival of large objects is limited by gravitational scattering during close plan- 
etary encounters and by collisions as well , if they spend sufficient time in the asteroid 
belt. The observed radiation-exposure ages of chondrites are shown to be consistent with 
this model . 


T he distribution of the masses of meteoroids 
is governed by several processes (Whipple, 
1967). Large numbers of new objects arc injected 
into the solar system by comets. Many small 
objects are removed by the Poynting-Robertson 
effect (Robertson, 1937; Wyatt and Whipple, 
1950) ; particles are destroyed by interparticle 
collisions and their shattering into fragments 
creates new particles (Whipple, 1967 ; Dohnanyi, 
1967). The influence of these collisions on the 
distribution of meteoroid masses has recently been 
discussed by Dohnanyi (1970; to be referred to as 
D-I in this paper) , who showed that the dis- 
tribution of meteoroids in the photographic range 
and of fainter ones is not likely to be stationary 
unless many of the particles, destroyed by colli- 
sions, are replaced by new ones given off by 


comets; the influence of radiation pressure on the 
size distribution of such cometary debris was also 
stressed (D-I). 

The orbital elements of meteoroids undergo 
frequent and random changes caused by planetary 
perturbations (Opik, 1951; Arnold, 1965). The 
influence of this process on the radiation closure 
age distribution of meteorites has recently been 
discussed by W 7 etherill (1967) and Wethcrill 
and Williams (1968). This age distribution was 
found to be sensitive to the survival times of 
meteorite-producing objects with respect to 
catastrophic collisions and to the rate at which 
these objects can “diffuse” through the solar 
system as a result of random gravitational per- 
turbations caused by close planetary encounters 
(Opik, 1966). 
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In this paper we shall discuss some aspects of 
the current evolution of the mass distribution of 
meteoroids; orbits will only be considered to the 
extent that they may influence the mass dis- 
tribution. This will be shown to lead to a self- 
consistent description of the dominating processes 
controlling the mass distribution of meteoroids 
with masses ranging from micrometeoroids to 
meteorite-producing objects having masses of tens 
of tons. 


OBSERVATIONAL EVIDENCE 

This section is a discussion of the observed dis- 
tribution of sporadic meteoroids ranging in size 
from micrometeoroids to large objects. 

The flux n(m) dm of meteoroids having a mass 
in the range m to dm incident on a unit area per 
2 Jr sr per unit time will be taken as: 


n (m) dm = am~ a dm ( 1 ) 

where a and a are constants in different mass 
ranges; a is known as the population index. 

Figure 1 is a plot of the cumulative flux N(m) 
of meteoroids into Earth’s atmosphere per meter 2 
per sec per 2ir steradians having a mass of m kg 
or greater. 



n(M) dM 


( 2 ) 


where M* is the mass of the largest object included 
among meteoroids. Near the small mass limit of 
the distribution, I used the results of the Pioneer 
8 and 9 data obtained by Berg and Gerloff (1970) , 
multiplied by two to correct for the Earth’s 
focusing effect. These authors have found an 
indication of a “cutoff” in the population of 
meteoroids, at a mass of about 5X10“ 15 kg, which 
sets the effective upper limit to the flux of pene- 
trating particles, as seen in figure 1. 

The points labelled Explorer XXIII and Pegasus 
are the influx rates measured by these satellites, 
and are based on calibrations by Naumann 
(1968) and Naumann et al. (1969). The penetra- 
tion sensors aboard Explorer XXIII and Pegasus 
were calibrated in the laboratory by firing particles 
at meteoric velocities into sensors similar to those 
actually flown. Since many of these particles w’ere 
accelerated gas dynamically, a fraction of their 



Figure 1— Cumulative flux (m" 5 s~ l /2ir sr) of meteoroids 
into Earth’s atmosphere having a mass of m (kg) or 
greater. 


masses may have ablated during the acceleration, 
so that the indicated flux is likely an upper limit. 
Some micrometeoroids may, however, be fluffy 
and of low density (Soberman, 1971) and would 
be less penetrating than were the laboratory 
particles of equal mass, so that the indicated flux 
may also be a lower limit. The nominal flux is, 
however, in agreement with the penetration flux 
measured by the Ariel II satellite (Jennison et ah, 
1967) . Since it is difficult to estimate precisely the 
uncertainties involved, an order of magnitude 
approximation may be the most accurate estimate 
that can be attained at the present time. 

Data from visual, radar and photographic 
observations as well as zodiacal light studies have 
been considered before (Whipple, 1967; Dohnanyi, 
1965). A best estimate to fit these data has been 
obtained in D-I; it has the form of equation (1) 
with 

a =1.5, for micrometeoroids, m<10~ 10 kg 



CURRENT EVOLUTION OF METEOROIDS 


365 


a= 1 % t for larger meteoroids but smaller 

than meteorite producing objects, 

10“ 10 kg<m< 1 kg (3) 

This simple model gives a good fit to the data, 
as can be seen from the appropriate portions of the 
curve in figure 1. Numerical flux values are given 
in table 1. 

Mass distributions obtained from satellite 
microphone measurements have been discussed by 
Kerridge (1970) and McDonnell (1971) ; because 
of calibration difficulties, many of these data are 
difficult to interpret. In those cases where calibra- 
tion difficulties have likely been overcome, the 
results are comparable to the penetration data 
(Kerridge, 1970), 

Mass fluxes of micrometeoroids estimated from 
particle collection experiments on board rockets 
and satellites are subject to uncertainties arising 
from contamination and identification difficulties 
(see Fechtig et al., 1968; Dohnanyi, 1971a, for an 
annotated bibliography) and wall not be employed 
in this study. Many of these particles are com- 
parable to or smaller than the wavelength of 
light; a discussion of the interaction of such small 
particles with the solar radiation field requires a 
discussion of interference effects and is beyond the 
scope of this paper. 

The influx rate of meteorite-producing objects 
has been estimated variously by Brown (1960), 
Hawkins (1963), and Opik (1958) ; their estimates 
are plotted in figure 1. These data are based on 
the mass distribution of recovered meteorites and 
their estimated rates of fall. Hartmann's (1965) 
estimate of the flux rate of large crater-producing 
objects is also indicated. 

Precise photographic observation of fireballs 
from the Prairie Network Project led McCrosky 
(1968b) to obtain as the cumulative flux (in 
MKS units) for these objects: 

N(m) = 10" 16 - 96 m -0 ' 62 1 <m< 10 4 kg (4) 

As can be seen from figure 1, this flux is about 
an order of magnitude higher than that of 
Hawkins' (1963) stones and about two orders of 
magnitude higher than the other earlier estimates. 
Extrapolation (dashed line in fig. 1) of the Prairie 
Network data leads to an even greater difference. 
Uncertainties in the photometric masses of these 
objects are not believed to span this discrepancy 


Table 1 —Differential Flux of Meteoroids n(m) =am a for 
Different Mass Ranges' 1 


Mass range a, m 3 kg* l /2v sr 


m<5Xl0 -ls kg <3/2 

5 X 10“ 15 kg <?7i < 10 -l ° kg 1.4X10 -11 3/2 

10 -10 kg <m <1 kg 3 XlO -18 13/6 


& Bee equation (1) in text. 


(McCrosky and Ceplccha, 1970); low density, 
fragile objects whose fragments do not survive 
atmospheric entry are believed responsible for 
the higher flux of fireballs than had been estimated 
earlier for meteorite producing objects. 


INFLUENCE OF COLLISIONS 

Meteoroids frequently undergo mutual colli- 
sions. Since these collisions are inelastic, the target 
particles may either lose a small portion of their 
mass (erosive collisions) or be completely broken 
up (catastrophic collisions). The net result is a 
change in the meteoroid distribution. 

The equation that expresses the dependence of 
the population on collisions can be written as 




dt 


dt 


+ 


dfjm, t) 
dt 


dm (catastrophic collisions 


5/(m, 0,i 

+ dm (creation by fragmentation 

dt 


(5) 


where f(mt) dm is the particle number density 
function, i.e., the number of particles per unit 
volume of space in the mass range m to m+dm. 
The number density /(m, t) dm is perpetually 
altered by erosive and catastrophic collisions and 
the creation of fragments in the mass range m to 
m+dm by the crushing of larger objects during 
inelastic collisions. 

It has been shown (Dohnanyi, 1969 and D-I) 
that, for a distribution with a population index 
a— 1% the contribution of particle creation, 
expressed by the last term in equation (5) , is 
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minor compared with the other two processes. 
The influence of collisions on the distribution of 
small particles is then dominated by catastrophic 
collisions, while erosion dominates the dis- 
tribution of large particles (D-I). 

The influence of catastrophic collisions on the 
particle population is readily calculated. We take 
for the number density of small objects in the 
mass range m to ra+dm at 1 AU from the Sun 

x 14 , 

/( w ) = 2 ^ n (™) ( 6 ) 

Here V m is the average Earth entry velocity and 
n{m) is given in table 1; the factor Y 2 is included 
to correct for the Earth’s gravitational focusing 
effect, the factor 4 results from averaging the 
velocity distribution over all directions and n(m) 
is given by table 1 . 

Using earlier results from photographic meteors 
(Dohnanyi, 1966) we take F« = 20.6 km/s and an 
average collision velocity for sporadic meteoroids 
equal to their mean geocentric velocity of 17.3 
km/s. In D-I, the largest meteoroid mass, r', 
that is catastrophically disrupted by impact of a 
projectile meteoroid of unit mass, was estimated 
(Moore and Gault, 1965; Moore and Robertson, 
1966). The result is r'«7.5X10 4 for basalt 
particles and about an order of magnitude smaller 
for pumice particles. 

Assuming a steady state distribution, one can 
readily calculate how many particles in this 
environment will survive catastrophic collisions 
after a time t . The result is shown in figure 2, which 
is a plot, for basalt-like particles, of the number 
density of particles that survive disruption after 
various time intervals, as indicated. 

It is readily seen, from figure 2, that the heaviest 
toll is taken from particles in the faint radiometeor 
range; less than 0.1 percent of these particles 
survive disruptive collisions in the mass range 
10 -9 kg<m<10 -5 kg during a time interval of 
about 10 5 yr. Larger particles survive longer 
because the number of projectile masses lethal to 
larger objects decreases; micrometeoroids, on the 
other hand, have a longer survival time for dis- 
ruptive collisions because the distribution of small 
objects that disrupt them tapers off. 

The true number of surviving micrometeoroids 
is smaller than indicated in figure 2, since the 
Poynting-Robcrtson effect will cause many of 



Figure 2 —Number density, /(m, t)dm, per meter 3 of 
meteoroids having masses in the range m to m+dm 
(kg) is plotted for different times /, as indicated, for 
particles surviving disruptive collisions in a stationary 
population /(m, o). 

these particles to spiral in toward the Sun. This is 
indicated as a function of particle mass in figure 3, 
which is a plot of meteoroid survival times with 
respect to removal and destruction processes. 

The survival time with respect to catastrophic 
collisions is made up of two main contributions: 

l/r=l/r cc +l/r, (7) 

where 

/ Afao 

f(M) dM (8) 

771 /T' 

and 
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1/T,(m) =fc [ M * M m f (M) dM (9) 

r cc (m) is the survival time of a particle with 
respect to fragmentation by a comparatively 
small projectile particle, and r s (m) is the survival 
time of a particle with respect to being swept 
up by larger objects. The Quantity K is defined in 
(D-I) as 

K=(3 7r i/2/4p)2/3{7) (10) 

where p is the material density of the particles and 
{V) is the average encounter speed. 

Figure 3 also shows a plot of the particle lifetime 
tpr with respect to the Poynting-Robertson effect 
for a particle of material density of 3.5 g cm- 3 
like that of basalt, and a material density of 
% g cm -3 which resembles pumice. 

The time for a particle to erode to M of its 
radius rs is plotted for a meteoroid with a density 
3.5 g cm -3 resembling basalt in composition for 
two linear erosion rates: 100 A yr~ l , which is 
an upper limit obtained earlier by Whipple 
(1967), and lA yr _I , a recent estimate based 
on cosmic ray track densities in glass removed 
from Surveyor 3 spacecraft, as w r ell as in some 
lunar samples (Fleischer et ah, 1971). 


The line labelled “gravitational lifetime” is the 
effective survival time of a particle in an Earth- 
crossing orbit with respect to being swept up by 
the Earth (Opik, 1966) . The actual time a particle 
is expected to spend in Earth-crossing orbit is 
considerably shorter (Wctherill, 1968); multiple 
near encounters give rise to a random wnlk process 
causing the particle to “diffuse” out of the region 
where its motion may be perturbed by the gravita- 
tional field of the Earth. 

Comparison of survival times with respect to 
the various dynamic processes plotted in figure 3 
indicates that the Poynting-Robertson effect 
determines the survival times for micrometeoroids. 
Also, the survival times for disruptive collisions of 
micrometcoroids with masses m < 10 -12 kg is 
dominated by the sw r eeping-up action of larger 
particles. Erosion is negligible for small objects. 
Catastrophic collisions with relatively small 
projectile particles dominate the survival times of 
particles in a mass range of 10" 10 kgCmClO 3 kg; 
larger objects will be dispersed by random w r alk 
from the vicinity of Earth's orbit long before they 
may be destroyed there by collisions. 

Survival times with respect to catastrophic 
collisions of asteroids in the asteroid belt 



Figure 3 . — Survival times as a function of mass of stray objects with respect to different loss 

processes, as indicated. 
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(Dohnanyi, 1969) and the extrapolated values of 
these survival times for small masses are indicated 
in figure 3 for the sake of comparison. 

INFLUENCE OF RADIATION PRESSURE 
ON THE DISTRIBUTION OF 
METEOROID MASSES 

A certain amount of radiation pressure, in a 
direction away from the Sun, is exerted by the 
Sun's light on all objects in the solar system 
(Robertson, 1937; van de Hulst, 1962). The result 
is a decrease, with decreasing particle radius, in 
the central force attracting the particle to the Sun. 
For particle radii smaller than a certain critical 
value, the electromagnetic force exceeds the 
gravitational force and the particle is blown out of 
the solar system. This critical radius is a function 
of the material density of the particle and its 
optical properties. 

Consequently, the size distributions of particles 
injected by a comet into orbits about the Sun with 
an initial angular momentum equal to that of the 
comet will have a “cutoff” below a certain critical 
size due to radiation pressure. This problem was 
discussed in D-I (also Harwit, 1963), and the 


results are summarized in figure 4. This figure is a 
plot of the critical particle mass just blown away 
by radiation pressure when ejected at the peri- 
helion of some parent comet having an eccentricity 
e, as indicated. The particles arc assumed opaque 
and particle material densities of I, 2, and 
3.5 X 10 s kg/m 3 have been considered. The eccen- 
tricities of some major showers are also indicated. 
One may assume that each comet that presumably 
gave rise to a major shower had, at the time it 
created the shower, an eccentricity similar to that 
of the shower. The intersection of each horizontal 
line, representing the eccentricity of the shower, 
with any line labelled by a density gives the 
smallest particle mass with the given density that 
can be present in the shower, according to geo- 
metric optics. 

It is readily seen, from figure 4, that none of the 
major showers would initially contain particles 
smaller than about 10~ 9 to 10~ 12 kg having a 
density of 3^ X 10 3 kg/ m 3 . The radiation pressure 
cutoff on the shower masses occurs in the range 
of ID” 11 to 10~ 13 kg for meteoroids having a density 
similar to basalt (3.5X10 3 kg/m 3 ). For particles 
having a mass of about 10 -15 kg or smaller, inter- 
ference effects come into play and a discussion of 



Figure 4. Radiation pressure limit for particles of different material densities p ejected by some 
parent object at perihelion; e is the eccentricity of the parent object. 
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the stability of those very small particles is beyond 
the scope of this paper. 

Similar considerations apply to fragments 
ejected from a parent body after an inelastic 
collision with another, smaller, meteoroid. In 
order to see, in detail, the effect under discussion, 
wu consider the eccentricity E and semimajor 
axis A of a fragment ejected from a parent body 
with eccentricity e and semimajor axis a. Following 
D-I we denote by v the ratio of the total central 
force on the fragment particle, i.e., the force of 
gravity reduced by radiation pressure, to the 
gravitational force acting on it. Since typical 
ejection velocities are small compared with the 
velocity of parent objects (Gault et ah, 1963; 
Dohnanyi, 1971b) we take the heliocentric velocity 
of the ejecta equal to that of the parent objects 
immediately after impact. Assuming conservation 
of angular momentum, it can then be shown that 

A = ra/[l — 2a(l — p)/ r 3 (H) 

where r is the radial distance from the Sun to the 
point of collision. Since p< 1, by definition, it can 
be shown from equation (11) that the semimajor 
axis, A, of the ejecta is always greater than the 
semimajor axis of the parent object. 

It is readily seen from equation (11) that if 

a/r>2/(l — v) (12) 

the ejecta wdll be expelled from the solar system 
by radiation pressure. 

According to equation (12), elimination from 
the solar system by radiation pressure is favored 
for collisions near perihelion (i.e., small a/r) and 
since the particle number density as well as the 
encounter velocity increases rapidly toward the 
Sun (D-I), we may expect most collisions to occur 
near the perihelion of the colliding objects and 
hence figure 4 may also be applied to collisions, as 
a rough approximation.* It may then be concluded 
that the contribution to the micrometeoroid 
population of fragments created during collisions 
by larger objects is much reduced by the elimina- 
tion of these ejecta by radiation pressure. 


* An exception occurs for objects having a node through 
the asteroid belt; depending on their orbital elements, it 
may then be physically possible for such objects to experi- 
ence a more severe collisional environment in the asteroid 
belt than at perihelion. 


The smaller the density of the fragment, the 
stronger is the effect of radiation pressure under 
discussion; this wuuld seem to create a natural 
selection favoring the elimination of fluffy 
particles. 

It has been show T n in D-I that a population of 
objects with a mass distribution having a popula- 
tion index of is not stable; because of collisions, 
particles in any mass range are destroyed faster 
than they can be replaced by the creation of frag- 
ments of the same mass range. Hence, if the present 
distribution of meteoroids has reached an approxi- 
mately steady state configuration, it is necessary 
to have a source of meteoroids replenishing the 
particles destroyed by collisions but not replaced 
by fragments. If meteor show 7 ers are indeed the 
required source which replenishes the population 
of sporadic meteoroids, and if the know T n major 
show T ers and comets are representative of the 
source of meteoroids, one wmild expect a drop in 
the population of sporadic _ meteoroids wdth 
masses smaller than 10" 9 to 10 -11 kg, depending 
on their density. This conclusion appears to be 
borne out by observation, as is indicated in the 
changing slope of the sporadic meteoroid mass 
spectrum around 10" 10 kg. 

RADIATION AGES OF METEORITES 

The survival time of large meteoroids with 
respect to catastrophic collisions increases mono- 
tonically with mass (fig. 3) ; the survival time of 
sufficiently large objects, therefore, will be 
dominated by the dispersive effects of gravita- 
tional perturbations. Taking the mean dispersal 
time of an object to be shorter than 1(F million 
years, w r hich is the lifetime of objects with respect 
to being swept up by the Earth, w 7 e see from figure 
3 that the survival time of objects with masses of 
hundreds of kg wall be limited by r<j. This, how- 
ever, does not necessarily hold for objects whose 
orbits cross the asteroid belt; depending on the 
relative time spent in the asteroid belt as w r ell as 
on the average encounter velocity with asteroids, 
the survival time of these objects may still be 
collision dominated. 

Figure 5 is a plot of the observed rate of fall of 
chondritic meteorites, as given by Wetherill 
(1969) in different cosmic ray exposure age ranges. 
It can be seen, in figure 5, that few objects have 



370 


EVOLUTIONARY AND PHYSICAL PROPERTIES OF METEOROIDS 



EXPOSURE AGES (MILLION YEARS) 

Figure 5. — Observed cosmic-ray exposure ages of chon- 
dritic meteorites (Wetherill, I960). 


ages smaller than a million years. Under steady 
state conditions it appears, therefore, that most of 
these objects originated in a region of the solar 
system far away from Earth’s orbit, and spent 
some time on their way to the Earth’s vicinity. 
Consequently, very few newly created meteorites 
are present. 

We therefore consider a simplified model for 
meteorite producing objects. We shall assume that 
the meteorite producing objects have orbits (cf., 
Wetherill and Williams, 1968; Wetherill, 1969; 
and Anders, 1971, for a more general discussion) 
with nodes that do not intersect Earth’s orbit, 
but do intersect the asteroid belt. This latter 
assumption is reasonable if we assume that the 
orbits of these objects resemble somewhat those 
of the Prairie Network fireballs (McCrosky, 
1968a). Using typical orbital elements for these 
objects, we have 

4 AU 

e~0.7 (13) 

10 ° 

where q is the aphelion distance, e is the eccen- 
tricity of the orbit and i is the inclination. It is 
readily shown that such an object spends about 
75 percent of its time at a solar distance between 
2.2 AU and aphelion. If all this time is spent 
within the asteroid belt, then use of a steady state 
asteroidal mass distribution (Dohnanyi, 1969) 
together with a mean encounter velocity with 
asteroids of about 11 km s -1 (as suggested by the 
orbital elements of equation (13)) leads, for an 


object with a mass of 10 4 kg, to 

r cc ~l million years (14) 

where r cc is the survival time of the object with 
respect to catastrophic collisions with asteroids. 
We chose a mass of 10 4 kg as a typical meteorite 
parent object; much smaller objects arc likely to 
produce few, if any, meteorites after atmospheric 
entry and larger objects are comparatively so 
scarce that their contribution to the production of 
meteorites will be neglected here. The value of 
1 million yr for t cc is an underestimate, since our 
test object is not likely to spend all of its time 
near the center of the asteroid belt; indeed, having 
an aphelion distance of 4 AU, it will spend about 
40 percent of its time beyond 3.5 AU from the 
Sun. Thus, taking 3.5 AU as the outer bound of 
the asteroid belt, 

T cc ~2*10 6 yr (15) 

appears a more reasonable estimate. 

The cumulative number of 10 4 kg objects, 
h(T), having an age of T million years or longer 
is then 

h{T) =hoe~ TlTe * (16) 

where ho is the total number of our meteorite 
parent objects with masses of 10 4 kg. It is assumed, 
in equation (16) , that these objects are in a steady 
state distribution and ho and h(T) are therefore 
independent of time. 

We now assume that these objects are occa- 
sionally perturbed into Earth crossing orbits with 
node(s) oscillating around 1 AU so that collision 
with the Earth becomes possible and the objects 
may then be recovered as meteorites. In order for 
this to happen, one or both of the nodes must be 
at a distance of 1 AU from the Sun which means 
that, for non-zero inclinations, the aphelion of the 
orbit will be some distance above the ecliptic and 
thus above the central region of the asteroid belt. 
The collision al survival time will be correspond- 
ingly lengthened for objects in such Earth- 
crossing orbits. As a first approximation, we 
assume, therefore, that these Earth-crossing 
objects will have survival times determined by 
the random gravitational scattering in close 
encounters with the Earth, i.e., the collisional 
survival times of these Earth-crossing orbits are 
relatively long. 
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We may then write 

dn/(t) _ n e '(t) h(t) ^ 

dt ~ y x 

where n/ (t) is the number density per unit volume 
of 10-ton objects in Earth-crossing orbit at a time 
tj and h(t) is the number density per unit volume 
of 10-ton objects in previous, non-Earth-erossing 
orbits. Quantity 1/y is the a priori probability, 
per unit time, for any Earth-crossing object 
to be scattered out of Earth-crossing orbit by a 
close encounter with the Earth and l/x is the 
a priori probability per unit time for a non-Earth- 
crossing meteorite parent object to be scattered 
into an Earth-crossing orbit. 

We assume steady state conditions: 

dn e '(t)/dt=0 (18) 

and hence 

rie=(y/x)h (19) 

as one would expect. 

To calculate the distribution of the ages T of 
these objects, we replace t by T in Equation (17) 
and get 

dn e (T)/dT=>-n e {T)/y+h{T)/x (20) 

where [dn e {T)/dT’] dT is the number density per 
unit volume of 10-ton objects with ages in the 
range of T to T-\-dT, and n*(T) and h(T ) arc 
the number of Earth-crossing and non Earth- 
crossing objects, respectively, per unit volume 
with ages greater than T* Because of the imposi- 
tion of steady state conditions, ?i e and h are given 
by equation (19), for T = 0, i.e., the number of 
objects for all ages is given by equation (19). 
Using equations (16), (19) and (20) we get: 

^4? -«.(0) (*-*•/'— «T«*) /(tw-y), r cc ^y 
dT 

( 21 ) 

and 

^4r^=n.(0 )T{<r T ^)/y\r cc ^y (22) 

61 

The expression for p».(r)/dT] dT determined 
by equation (21) or (22), gives the number of 
meteorite producing objects having an age of T to 
T+dT and is therefore proportional to the 
number of meteorites in that cosmic ray exposure- 
age range. It is readily seen, from equations (21) 


and (22), that the number of "young” meteorites 
having ages much smaller than y and r cc is zero. 
This happens because a finite time is required for 
an object, which has just been created in a non- 
Earth-crossing orbit, to “find its way” to an 
Earth-crossing orbit. 

A comparison of equation (21) with the observed 
cosmic ray exposure ages of chondritic meteorites 
based on a compilation by Wetherill (1969) is 
plotted in figure 6. The data published by Wetherill 
(1969) appear in the form of histograms (fig. 5) 
representing the rates of fall for meteorites in 
various ranges of cosmic ray exposure ages. In 
figure 6, WcthcrilFs data are reduced to represent 
the differential distribution d?i e (T)/dT of cosmic 
ray exposure ages per unit exposure-age range. 
Taking 

tte(0) =68.95 (23) 

which is the total rate of fall for all exposure ages 
given by Wetherill (1969), dn e (T)/dT from 
equation (21) is then plotted in figure 6 for several 
combinations of the survival times r cc and y. 

It can be seen from equations (21) and (22) 
that d?i e (T)/dT is symmetric with respect to t cc 
and y. This means that if either r ec or y is very 
short, there will be a peak in the distribution of 
meteorites with very short exposure ages. This 
happens because if r cc is short, young objects are 
favored and if y is short, meteorites do not have 
enough time to “grow old” in Earth-crossing 
orbits. If, on the other hand, either r cc or y is 
long, the distribution has a long exponential tail 
containing many old meteorites. This happens 
because relatively many objects escape cata- 
strophic collisions if r cc is long, which results in a 
corresponding abundance of old objects. If y is 
long, then meteorites in Earth crossing orbits stay 
around long enough for many of them to age 
there. If r cc is comparable to y } then most objects 
have an age of that same order of magnitude, 
i.e., a narrowly peaked distribution of radiation 
ages. 

It can be seen, from figure 5, that the best 
estimate of 2X10 6 yr for r cc (cq. (15)) combined 
with a gravitational survival time in Earth- 
crossing orbit of about 2X10 6 yr provides a 
reasonable fit to the observational data. 

We now consider the frequency of meteorite 
falls. According to the results of the Prairie 
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Figure 6. — Observed and theoretical distribution of cosmic ray exposure ages per unit- exposure- 
age range of chondritic meteorites for different survival times t cc and y discussed in the text. 


Network Project (McCrosky, 19G8b), the flux of 
objects per unit mass range having a mass of 10 
tons that enter the Earth’s atmosphere is 
2.25 XlO -24 kg -1 m -2 s ~ l . Using equation (6) and a 
mean Earth entry velocity of 20 km s _l , we obtain 
a number density of 2.25 X 10 28 particles/ (kg m 3 ) . 

We now assume that the original meteorite 
parent objects in non-Earth-crossing orbits are 
objects in Mars-crossing orbits. For the number 
density per m 3 kg of Mars-crossing objects having 
a mass of 10 4 kg, we employ the results of Kessler 
(1970), who estimated the number density of 
large objects at a distance of about 1.5 AU from 
the Sun to be about 1 percent of the peak number 
density of the asteroidal belt. Assuming the latter 
to be about double the average number density 
(Narin, I960) and using the results of Dohnanyi 
(1969), we then get a value of 2.22 X 10~ 28 for the 
number density h per (kg m 3 ) of 10 4 kg objects at 
a distance of 1.5 AU from the Sun (i.e., Mars- 
crossing orbits). This number is, however, an 
underestimate, since these Mars-crossing objects 
spend most of their time away from perihelion and 
near aphelion. A best estimate for the number 
density of Mars crossing meteorite parent objects 
is perhaps an order of magnitude greater. 

Using equation (19), we can now calculate the 


survival time x; the result is 

x = y(h/n e ) « 10X2(2.22 X 10- 28 /2.25 X 10~ 28 ) 

«2Xl0 7 yr (24) 

which is surprisingly short. On the other hand, if 
we used Hawkins’ (1963) estimate for the flux of 
stone meteorite-producing objects, then the num- 
ber density n e of 10 4 kg objects near Earth would 
be 6.4XlO“ 30 /(kg m 3 ) and we would have 

x = 2X 10(2.2 X 10 _2S /6.4 X 10“ 30 ) 10 6 yr 
= «7Xl0 8 yr (25) 

It therefore appears that a diffusion time x 
ranging from as short a value as 2 X 10 7 yr up to 
about 6 XlO 9 yr can explain the data, depending 
on the flux law employed, or more specifically, 
depending on the extent meteorite producing 
bodies are assumed to contribute to the Prairie 
Network flux of fireballs. 

DISCUSSION AND CONCLUSION 

The considerations of the section on radiation 
pressure suggest the difficulty of creating micro- 
meteoroids in substantial quantities, either by 
cometary emission or by the production of second- 
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ary ejecta during collisions: Radiation pressure is 
likely to expel most of the small objects shortly 
after their creation, while those that initially 
survive radiation pressure vail spiral into the 
Sun as a result of the Poynting-Robertson effect. 
While a detailed mathematical formulation of this 
problem has not been developed, the “bending 
over” of the flux curve for micrometcoroids (fig. 1) 
is probably caused by this process. 

Meteoroids in the radio and photographic 
meteor range have survival times limited by 
catastrophic collisions. Fragments from the dis- 
ruption of relatively large objects are not pro- 
duced in sufficient quantities to maintain the 
population of these meteoroids in steady state 
(D-I) , hence it is necessary that a steady supply 
of cometary debris be available to replenish the 
destroyed particles. 

Large objects, of the meteorite producing class, 
have survival times limited by collisions as well as 
gravitational dispersal. If the latter were un- 
important, most meteorites would have a very 
young radiation exposure age corresponding to 
the curve in figure 6, labelled (r co y) = (2, .01) 
m.y.; also, if collisions were very rare and the 
survival time were dominated by gravitational 


dispersal alone, a similar exponential type curve 
would be obtained. The data support neither of 
these distributions but suggest that the gravita- 
tional dispersal time is comparable to the eolli- 
sional survival time of these objects prior to their 
having been scattered into Earth-crossing orbits. 
This is consistent with a model in which the 
meteorite parent objects have Mars-crossing 
orbits and aphelia in the asteroid belt or some- 
what beyond (Wetherill, 1968, 1969; Anders, 
1971). These objects are then scattered, during 
close encounters with Mars or perhaps Jupiter, 
into Earth-crossing orbits. This simple model is 
found to predict a distribution of meteorite 
radiation exposure ages consistent with observa- 
tion (Wetherill, 1969). While a precise identifica- 
tion of these original objects is still under dis- 
cussion (Wetherill, 1969; Anders, 1971), results 
derived from this simple statistical model may 
provide additional clues concerning their identity. 
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